VORTEX LIQUIDS AND THE GINZBURG-LANDAU EQUATION 

MATTHIAS KURZKE AND DANIEL SPIRN 



Abstract. We establish vortex dynamics for the time-dependent Ginzburg-Landau 
equation for asymptotically large numbers of vortices for the problem without a 
gauge field and either Dirichlet or Neumann boundary conditions. As our main 
tool, we establish quantitative bounds on several fundamental quantities, including 
the kinetic energy, that lead to explicit convergence rates. For dilute vortex liquids 
we prove that sequences of solutions converge to the hydrodynamic limit. 



Contents 



1. 


Introduction 


m 


2. 


Notation, assumptions, and some identities 


M 


3. 


Measuring true vortex positions 


M 


4. 


Initialization of Gronwall argument 


M 


5. 


Quantitative rate of second-order gamma convergence 


m 


6. 


Quantitive bounds on the kinetic and excess energy 


ESI 


7. 


Quantitative bounds on the supercurrent 


m 


8. 


Completion of Gronwall argument 


m 


9. 


Estimates on the Neumann function and the renormalized energy 


m 


10. 


Hydrodynamic limit 


m 


References 


m 



1. Introduction 
Let n : i7 — )• C satisfy the scaled Ginzburg-Landau equation 

(1.1) -^dtu = Au+ \u(l-\u\'^) 

\loge\ ^ ' 

with either Dirichlet boundary conditions 

(1.2) u = e^"^+^'^- on OQ 

with ipi, E C^, Jq^ dr^i, = 0, so deg(n; d^) = n, or Neumann boundary conditions 

(1.3) djyu = on dn. 
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We take to be a smooth, simply connected domain containing the origin. Equa- 
tion (1.1) models the dynamic behavior of superconductors when the electromagnetic 
field potential is absent. When a gauge field is present, the corresponding Gorkov- 
Eliashberg equations 

2 „, , ^ „. fi l„.l2\ 



(1.4) ■ ^" ' £2 

E = — curl curl A + jAiu), 

where = dt + i^, E = dtA + V<I>, and jA{u) = {iu, Vau), provide a more complete 
model of superconductivity. 

In order to describe the behavior of solutions of (1.1) with small e we define some 
fundamental quantities including: 

(1.5) energy density ee(n) = - iVnl^ + — r (l — luP)^ 

(1.6) supercurrent j{u) = {iu,Vu) 

(1.7) vorticity/Jacobian J(n) = det Vu =- curl j(u). 

Here (•, •) denotes the real scalar product of two complex numbers, so (o, b) = ^{ab + 
ah) for a,6 G C. Solutions to equation (1.1) diffuse the Ginzburg-Landau energy 

(1.8) Eeiu) = I e,iu) 

Jn 

via the identity 

(1.9) E,{u{t)) + f [f^ = E,{u{0)). 

Jo Jn |loge| 

1.1. Vortex dynamics and vortex liquids. A prominent feature of type II super- 
conductivity is the presence of localized regions, called vortices, where superconduc- 
tivity vanishes. In particular there exist some points {oj}"^^ in Q where \u{aj, •){ 0. 
Furthermore, about each vortex the winding number of the phase is quantized; in par- 
ticular 

^ [ T.j{u)^deZ\{0}. 

In the vicinity of each vortex the Ginzburg-Landau energy Es (u) blows up at the rate 
TT I log e I +0(1). Bethuel-Brezis-Helein showed in [3] that minimizers of the Ginzburg- 



Landau energy (1.8) can be expanded further up to second order 

E,{u) = n{7r\log e\+-fo) + W (a) + 0{l), 
where 70 is a universal constant and 

W{a) = —IT log \ ai — aj\ + boundary effects 

is a renormalized energy and the winding number about each vortex is one. This 
renormalized energy is precisely the bounded domain version of the Kirchhoff-Onsager 
functional that arises in two dimensional incompressible Euler equations and other 
settings. The renormalized energy will be discussed in more detail in Section [2] and 
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Section |9j From back-of-the-envelope calculations one finds that J{u) is quantized 
and looks like a sum of integer-weighted delta functions; and so, for small e one finds 
that 

|loge| ^ 

in case when the winding number about each vortex equals one, and as e — t- 0, u 
limits to 



(1.10) = n 



where V'* is H^{0,). This is referred to as the canonical harmonic map when ip-i, is 
a harmonic function. This limiting behavior was established in many situations, see 
for example H |32l 1121 [M [20] . 

When dynamics (1.1) are turned on, these vortices move according to the gradient 
flow of the Kirchhoff-Onsager energy: 

(1.11) dj = --Va,W. 

vr 



The I log e I factor in front of (1.1) is the critical time scale on which vortices will move 
and can be thought of as the length of time it takes the unsealed time dependent 
Ginzburg-Landau equation to move an 0(|loge|) amount of energy an 0(1) distance. 
That vortices satisfy ( |1.11[ ) in the limit was the subject of a formal asymptotic study 
by E [12]. Later, arguments of Lin [30] and Jerrard-Soner [21] provided rigorous 
justification of the limit. Both [30] and j21j assume that the number of vortices is 
uniformly bounded as e — )• 0. The limit equation ( 1.11 ) is the gradient flow of W just 
as (1.1) is the (rescaled) gradient flow of the integrated energy density J^ee{u). The 
similarity in structure can also be seen by the energy dissipation identity 

(1.12) W{a{t))+ [ \h{s)f = W{a{0)). 



This structure was exploited to give a more abstract proof of the motion law by 
Sandier-Serfaty |l3] in their F-convergence of gradient flows framework. 

In recent years there have been significant advances in understanding the dynamics 
of a finite numbers of vortices by Bethuel-Orlandi-Smets [3] on and by Serfaty [48] 
on bounded domains. These results allow for much weaker initial conditions, handle 
collisions of plus/minus vortices, and describe the dynamical behavior of higher degree 
vortices. 

On the other hand, the behavior of the time dependent Ginzburg-Landau equations 
with asymptotically large numbers of vortices has received mostly formal treatment. 
The question of how large numbers of vortices behave in superconductors is important 
from both experimental and numerical perspectives. In the former, typical supercon- 
ductors contain many millions of vortices per sample [6l[T3] so the large vortex problem 
is a fundamental feature of high Tq superconducting devices. In the latter, point vor- 
tex methods provide a useful class of numerical algorithms for simulating challenging 
PDE's, like vortex sheets; hence, ( |1.11[ ) is a reasonable numerical approximation of 
the limiting mean field equation with vortex sheet initial data. 



4 



M. KURZKE AND D. SPIRN 



In |13j E looks at how the planar analogue of (1.11) behaves in a mean field sense 
as 7T, — )• oo. Defining the vortex density function a;„ = - Yl^=i ^aj(t)^ the author shows 
that the limiting density, oj = lim„_>oo '^m formally satisfies a weak PDE of the form 

dtoj + d\v{ijjv) = 

u = V(A-^)a; 

after rescaling time t. Subsequently, this ODE limit on was rigorously established 
by Lin-Zhang [34J. Furthermore, these authors show that this Euler-type equation 
has many regularizing features that will be discussed later. 

There are many similarities between this ODE limit problem and ODE limit prob- 
lem arising from the point vortex method for the Euler equations. In the latter case 
it was shown by Schochet [46] and Liu-Xin [36] that the vortex density function for 
Euler point vortices on M^, which follow the Kirchhoff law 

aj = --Vj^W{a), 

limit to a weak Delort solution to the incompressible Euler equations on M^. Due to 
the similarities of the problem, Lin-Zhang |3l] use the approach of [36j to prove the 



associated hydrodynamic limit of (1.1) on planar domains. 

The present work is the first to directly couple the Ginzburg-Landau equation to a 
mean field PDE. All previous works either prove a PDE to ODE limit for a bounded 
number of vortices or an ODE to mean field PDE limit. Our quantitative results 
enable us to take the diagonal limit in a rigorous way. 



Ginzburg-Landau equation 




Vortex motion law 


► 








n— >" 



IVIean field equation 



Figure 1. Limiting from the Ginzburg-Landau equation directly to 
the mean field equation 

In order to make the direct connection between the Ginzburg-Landau equation and 
the limiting mean field equation, it is necessary to establish two steps. The first of 



which entails a proof that (1.1) can accept asymptotically large numbers of vortices 
for long-enough times. The second step involves coupling these Ginzburg-Landau 



solutions to an appropriate hydrodynamic limit of (1.11) on bounded domains. 



1.2. Results. We now state our first theorem which establishes vortex dynamics for 
asymptotically large numbers of vortices. Though the initial conditions require well- 
preparedness and all degree plus-one vortices, such conditions can be softened a bit 
by combining with the work of [H |38] . In the following we let 

A<B \iA<CB 
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for some C that depends only on Vt. Next we define the excess energy 

D{a{t)) = Ee{u{t) - [n(7r|loge| +7o) + W^(a(t))] , 
which will be used to control the deviation of the vortex path from the path defined 



by the ODE \l.n\ . We also define 

Pi, = - min {min I aj(t) — ak{t) \ , mindist(a,(t), c?r2)} 

to be a measure of how close vortices can approach each other or the boundary on 
the set time scale. Finally, we introduce a weak topology related to the length of a 
minimal connection, see |5], 



sup 

l|V<^|liOO(o)<l 



This norm provides a good scale-invariant measure of the distance of J{u) and -S^i^ 
to a sum of delta functions. In particular if |a,- — bj\ < for j = 1, . . . ,n then 



|log£| 



j j 



Y^\aj-bj\. 



We can now state our first theorem which supplies a long time existence result for 
asymptotically large numbers of vortices in the dilute regime. 



Theorem 1.1. Suppose u solves (1.1) with either (1.2) or (1.3). Furthermore, let 

1 i_ — — — 

n < I log e| 200 and p^, > \loge\ loo and 

(1.13) D{a{0)) < llogepi 

n 

(1.14) \\J{u{0))-2nY,KM\\w-^'Hn)<C.\logen 

then for all < t < T we have 
„e.(«)(t) 



|loge| 



\\Jiu)it)-Y^ 



/ 



Qp^ (a) 



eei\u\) + 



1 



4 

\u\ 



\u\ 



^ |loge| 4 
_ 1 

^ |loge| " 

^ llogep^ 
< |loge|~5i5 



where 
(1.15) 



T = min{CVlog|log£|^,ro} 

n 
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where tq = info<t<T{Pa(t) < P*}} solve and C = C(r2). Here ^p{a) = 

From Lemma 14 of |23j one can easily construct maps ul{x;a) that satisfy the 
well-preparedness assumptions ( 1.13 )-( 1.14). In the case of a bounded number of 



vortices is well known that this hypothesis is not very important, since one can show 
that data will become well-prepared almost instantaneously due to strong convergence 
estimates, see [30^ [21} HI HH], and we have no reason to expect a different behavior 
here. 

Given the result above, we can prove that the sequence of solutions converge in a 
prescribed sense to the expected hydrodynamic limit. In this theorem we study only 



Dirichlet boundary conditions (1.2) since we need to have the vortex motion law hold 
for times of order 0(n~^), and in the Neumann case (1.3) vortices will migrate to the 
boundary too quickly. 

What type of equation do we formally expect the vortex density to satisfy in the 
limit? If we rescale time t = nt then the limiting vortex density function uj = 
limn-5.00 ^ Sj=i ^aj{t) on a bounded domain should satisfy a weak form of the system 

djU) + div{ujv) = 

v = V 



(1.16) 



'Ax') 



where A_^^ 



a 



(1.17) 



w arises through the Poisson problem 
Alt; = g vaVt 
d,,w = d-rO on dVL 



for = arctan(y/x). Note consistency requires ^^g = 27r due to the Neumann 



boundary condition. To motivate a notion of an interior weak solution of (1.16) we 



follow Lin-Zhang [34]. If a; is a smooth solution to (1.16), we multiply by x £ Co°(^) 
and integrate by parts. Then — // dtX'^ ~ 1 1 (^x^XVkdxjVj = 0, where we used 
divf in the interior of Q,. Finally, using diV2 — d2Vi = and integrating by parts 



UJ 



again yields (1.18) below. We note that this definition is similar to the one introduced 



in [HI [38] for weak solutions to the Euler equations except that the associated test 
functions are exchanged. 



Definition 1.2. We say uj is a generalized interior weak solution to (1.16) if for all 



'1-^2} + I [ 29,1 5,2X V1V2 = 
10 Jn 



(1.18) - / / ujdtx + 

Jo Jn Jo Jn 

where Vj{x) = d x^Aj^ uj = J dx^N{x,y)uj{y)dy where N{x,y) is the Neumann func- 
tion defined via (1.17). 



We can now state our main result which shows that we can solve (1.16)-(1.17) with 



vortex sheet initial data via a sequence of point vortex methods with appropriate 
data. 

Theorem 1.3. Assume that ujq ^ -M. Ci H^^{il.) satisfies loq > 0, f^uJo = 27r, and 
supp(wo) C {dist(x, d^l) > Co > 0} for some constant Cq Then there exists a sequence 
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of initial data u^^ (0) wi th n = | log | log | log Enl | h number of vortices that satisfies the 
hypotheses for Theorem ^l . 1^ such that ^^^^^p^^^y^ 



ujQ in AinH as e„ — )• 0. Such 



(1.2) for times up to T 

nt and letting uJs{t) 



initial data generates a sequence of solutions u^^{t) of (1.1) with boundary condition 

log| logenll^- 

ee("e(t)) 



Setting t 



then 



n 7r|loge| 

oje(t) uj(t) in MnH~^ 

where uj{i) is a generalized interior weak solution, defined above, to 

dfOJ + div (vu) = 



(1.19) 



V A 



for all t G [0, oo). Finally, v{t) G Lf^^{Q) for all time. Here f = w if 

Aw = f inVt 
dyW = drO on d^l 



(1.20) 



where 9 = arctan(|) and f = 2-k. 

1.3. Discussion. The issue of whether the weak solution satisfies the correct bound- 
ary condition is a deep and difficult question. Since vorticity can (and should) concen- 
trate on the boundary, it is difficult to acquire the necessary regularity to ensure the 
boundary conditions are achieved in the classical weak sense. Some recent progress 
has been made in [T7| by establishing boundary- coupled weak solutions of the two 
dimensional incompressible Euler equations in exterior domains. 

To make a fully consistent limit it would be interesting to study the question of 



uniqueness of the limiting mean field equation (1.19). In [34j the authors establish 
uniqueness for initial data in L°° with compact support for the planar problem. A 
similar study of regular solutions would be natural for ( 1.19| )-(1.20) too. 

From formal considerations of (1.19) the vortex density function satisfies dtto + 

,2 



V ■ Vlo = —UJ , so along the trajectory of the induced velocity one sees that the 
density function should decay like t^^. In fact for smooth initial data on Lin- 
Zhang [34] proved this fact, which implies that the vorticity spreads out quickly from 
a compact set. This behavior implies that we expect most vortices to be pushed to 
out the boundary in a similar fashion. This conforms to the picture presented in 
Sandier-Soret |45] for global minimizers of the functional E^{u) on bounded domains 
constrained to the boundary condition of the type u = e*"^ and n — t- oo. Sandier- 
Soret show that vortices accumulate close to the boundary of the domain as n grows 
asymptotically large. Taken together, we should view Theorem 1.3 as a mean field 
description of the vortex density for times in the mesoscale in the interior of the 
domain. 

The dilute density of the vortex liquid results from two issues. The first is that we 
use energy comparison and a Gronwall inequality to pin the vortex positions to the 
ODE (1.11). This results in an upper bound T < —\/\ log |loge| | in Theorem 
Integrating methods of [l8] and/or [4J should improve some of these bounds. 



1.1 

te^ 

second issue arises from the poor bounds on the intervortex distance for the ODE 
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( [l.lip . Better knowledge of how the ODE behaves should improve the vortex density 
allowed here. 



Although (1.1) provides a fertile ground to test the mathematics of the Gorkov- 



Eliashburg equations, the more physical problem entails looking at the hydrodynamic 



limit of (1.4). For the Gorkov-Eliashberg equations (1.4), corresponding proofs of the 
vortex motion law are due to the second author [51j for 0(1) fields and Sandier- 
Serfaty |43j for larger fields, following the formal asymptotic work of [IQ]. Formally, 
it was shown by Chapman-Rubinstein-Schatzman [7] that the hydrodynamic limit of 



the associated ODE arising from the vortex motion law of (1.4) converges to a weak 
solution of 

dti^ + div(a;u) = 

V = V (A-iy^ uj. 



A rigorous approach of the hydrodynamic limit of (1.4) will be studied in [29]. 

For these problems, there has been a lot of recent progress on the limiting equations 
for the vortex densities. Ambrosio-Serfaty [1] and Ambrosio-Mainini-Serfaty [2] study 
them as a metric gradient flow in the space of measures, with the Wasserstein distance 
as the natural metric. However, they do not obtain the convergence. Even when it 
becomes possible to carry out the program outlined in the survey of Serfaty |49j and 
to directly obtain the Wasserstein gradient flow studied in [H [2] from the Gorkov- 
Eliashberg equation by a F-convergence of gradient flows type result, we believe that 
our approach will still be useful. For one, it provides quantitative bounds that are 
useful in type II superconductors without going to the e — )• limit of "extreme" type 
II superconductivity. More importantly, as our approach does not rely on the gradient 
flow structure, it can be adapted to yield results for more general situations, such as 
the mixed flows studied in [26] and [39] for the ungauged problem and in [27] and [50] 
for the gauged problem. Such motion laws have physical importance, as they can be 
used to explain the sign-change in the Hall effect of type II superconductors, see [TT] . 
|24) . Similarly, we expect that our approach can be adapted also to the Hamiltonian 
Ginzburg-Landau wave system, where results for the PDF to ODE limit for finitely 
many vortices have been found in |18] and [31], and the ODE to mean field PDF limit 
has been studied in [351. 



1.4. Method. We finish the introduction with an outline of the arguments in the 
paper. The general scheme of the paper is to carefully deduce the vortex motion 
law for the time dependent Ginzburg-Landau equations by energetic methods. These 
methods, which in many ways are not as strong as the elliptic PDF approach of Serfaty 
[171 SH] or the parabolic PDF approach of Bethuel-Orlandi-Smets do provide a 
way to avoid using convergence properties too early in the proof. In particular a 



fundamental feature of the arguments of all previous vortex motion law proofs in ( 1.1 ) 
is to use a limiting kinetic energy lower bound. This bound proves to be difficult to 
carry to the problem of asymptotically large number of vortices. Furthermore, though 
our setting provides weaker control on the rate of convergence of u to since we do 
not take advantage of strong convergence; nonetheless, we are still able to establish 
rates of convergence and long time of existence at the finite e level. Furthermore, 
since we are less reliant on strong convergence, our approach is more applicable to 



VORTEX LIQUIDS 



9 



situations where parabolic/elliptic regularization fails, such as in the mixed Ginzburg- 
Landau equation |iog£| dfU+idtu = Au+-^u (l — |np) and the Ginzburg-Landau wave 
equation |iogg| dfu = Au + (l — Our approach follows the program of the 

second author and R. Jerrard |23| for the Gross-Pitaevsky equation idtu = Au + 
p-M (l — Surprisingly implementing this approach for is more challenging 

and requires several new estimates including an explicit estimate relating the kinetic 
energy of the PDE to the kinetic energy of the ODE in Section |6] that are not needed 
in the Gross-Pitaevsky problem. 

In Section [2] we provide some background on the renormalized energy and discuss 
the differential identities on \u\, j{u), and Ce^u). Finally, we end the section with a 
list of some asymptotic bounds on fundamental quantities. These asymptotic bounds 
will either be prescribed through the choice of initial data or arise as a consequence of 
the results in the paper. Though the asymptotic bounds seem somewhat arbitrarily 
chosen, they represent a class of data for which we can close the arguments. Though 
there is room for improvement in the choice of these asymptotic rates, the order of 
magnitude of each choice is constrained by our arguments. 

We then turn our attention in Section [3] to a discussion of how to pin the vortex 
positions in (1.1) to the ODE positions. At each time instant one can prove, via 
localization estimates discussed in Section jsj that is close to a sum of delta 

functions with the error controlled by the amount of the excess energy. However, 
despite knowing the location up to a bunch of point locations ^j, it is not known a 
priori how close the ^j's are to the a^-'s, the expected ODE locations. One could 
try to establish control on a quantity like X]j=i ~ ^ Gronwall inequality. 

Though this approach works in the e — )• limit, see for instance |M1 El ESI EI], it 
does not work well for the fixed e problem. In particular one does not know if the 
^o 's are Lipschitz at the fixed e level. As in [23] we then introduce a weaker measure 



V = I |'iog£| X(^ ~ aj){x — aj)\ for a suitable cutoff function x- This function has 
the same effective behavior but is differentiable. At this point we define a sequence 
of timescales on which r] is well behaved in different senses. On these timescales it 
is shown that r]{t) is close to ~ along with other canonical measures. 

Therefore, our game plan is to prove that rj{t) is small on these define timescales via 
a Gronwall inequality; and hence, pin the vortices to the a^'s. 

Section [4] begins the Gronwall argument by computing ^rj, and this computation 
yields several terms that can be directly estimates found in Sections [5]|6j We remark 
that this line of attack works only for some of the terms at arise from the time deriva- 
tive of rj. Several terms arise that can only be estimated by which creates havoc 
within the Gronwall argument. Rather we will later look at a Gronwall argument for 
a time- averaging of r/, and this proves to be sufficient when we combine this with a 
Lipschitz bound on r] proved at the end of this section. 

Section [5] contains several technical results, a few of which arise out [23], on the 
Ginzburg-Landau energy. The first result state that if J{u) is close to a sum of delta 
functions located at a = {oj}, then we can find an even better collection of vortex 
locations ^ = {^j} that localize J{u) up to a small error that depends on the excess 
energy. The second result in this section provides a rate of convergence of u to 
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u^, = n^(a) outside of the vortex cores with an exphcit errors that depends on the 
locahzation of the Jacobian to the a = {ctj} and the excess energy. These first two 
results are established in [23j based on refined Jacobian estimates and depend heavily 
on results from [22]. The third result is a new estimate on the localization of the 
Ginzburg-Landau energy density to the same set of delta functions 



|loge| 



n 



< 



1 



loge| 



[nD{a) + W{a)] 



VF-i'i(n) 

The estimate presented here is a refined (i.e. e-rate dependent) version of an estimate 
found in Colliander-Jerrard [8]. 

Since the localization and gamma stability error estimates depend explicitly on the 
excess energy, it is necessary to understand how the excess energy evolves in time. 



By the energy dissipation identities (1.9) and (1.12) we see that, 

D{a{t)) = D{am + f \a\' - f f 

Jo Jo 7r|loge| 

consequently, D{a{t)) can be controlled by well-preparedness of the initial data and 
a lower bound on the kinetic energy. The main result of Section [6| and hence on the 
excess energy, is Theorem 6.1 which provides the lower bound 

\dtu\^_ _ /■* 

'0 



I ■ |2 \ 

a > 



-error. 



p-k, T, and supo<t<j' This result pro- 



10 7r|loge| 

where error depends explicitly on e, n, 
vides a purely quantitative approach to the kinetic energy lower bounds that are 
found in [31^ [T8l H3j, each of which rely on compactness properties to get a lower 
bound. To establish this result we make quantitative the kinetic energy estimate of 
[HTj , who used the differential identity for the conservation law along with an equipar- 
tition of energy result. Here we make use of an optimal equipartitioning result in |28| 
that identifies how close the tensor Vu Vu is to the diagonal matrix \dju\'^Sjk {djk 
being the kroenecker delta function). Placing this equipartitioning result into the 
differential identity for the Ginzburg-Landau energy ee(n), applying a test function 

xix — aj{t))aj{t) ■ {x — aj{t)), x a smooth cutoff function, 

and integrating over yields the lower bound. 

Though we are able to estimate most terms arising in our differential identity for 
^r/ in Section |4] by 77 or by a small error depending on e, there are several terms 

that have the general form Q ~ •^(^*)) some explicit function Q. A naive 
bound on such terms yields an estimate of the type which, unfortunately, destroys 
the sensitive Gronwall argument. Instead we prove that these terms are small directly 



by showing 



j{ui,) is small on the time scale that we are considering. This is 



achieved by employing a Hodge decomposition ^^j^ — j(n*) = V/i -|- V-^ f2- We show 
/i is small by the localization estimat on the Jacobian and /2 is small via an involved 
time-averaging argument. This method was employed in [23j in the Gross-Pitaevsky 
equations via a time-averaging technique. In Section [7] wc prove an analogue of this 
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type of bound for (1.1) for the parabolic problem. Unfortunately, since we have a 
different differential identity for j{u) and weaker bounds on many intermediate results, 
we feel it necessary to prove this proposition in detail for the sake of completeness. 

Finally we complete the proof of the motion in Section |8j This entails showing 
that we can extend the time-scales introduced in Section [2] to a long enough time 
T claimed in Theorem |1.1[ In particular we show by a Gronwall argument that the 
time-averaged rj behaves well on the long time scale, so long as vortices do not get 
too close. 

At this point we turn our attention to the hydrodynamic limit. Our first priority is 
to build a suitably large class of initial data for which there are asymptotically large 
number of vortices n — t- oo as e — t- and for which this class of initial data satisfies the 
ODE for times on the order of at least 0{-). We note that for (1.1) with Neumann 
boundary conditions, vortices will not remain away from the boundary of the domain 
for long enough times. Therefore, we turn our attention to the Dirichlet problem in 
Section [9j By a careful consideration of the renormalized energy W{a) we show that 
vortices obey the same logarithmic repulsion at the boundary that they feel from each 
other. This implies a long enough time of existence of the ODE ( |1.11 ). Initial data 
u{0) with vortices satisfying the conditions prescribed for our class of ODE initial 
data can then survive on a long enough time scale for our hydrodynamic limit. 

In Section 10 we complete the proof Theorem 1.3 on the hydrodynamic limit. The 
first part of the section we proof an analogue of results of |46t [36 j for the Euler point 
vortex method on Mp and the gradient fiow version of [Mj on for bounded domains. 
By a careful expansion of the time dependent behavior of ^Yyj=i^aj{t) integrated 
against a test function with compact support. By using estimates on the Neumann 
function and implementing the strategy of j36t [M] we prove the convergence and the 
local energy bound. To complete the proof of the theorem, we show that nonnegative 
vortex sheet initial with compact support can approximated by a sequence of a sum 
of degree-one vortices that satisfy the conditions of our class of initial data. The novel 
approximation scheme leads to an asymptotically long time of existence with asymp- 
totically large numbers of vortices. Finally, we set a;„(t) = ^ "^iiogei^ ' which has similar 



behavior as the vortex density function ^ X]j=i ^aj(t) due to Theorem 1.1 Relating 
this density to the ODE hydrodynamic limit, then achieves the desired hydrodynamic 
limit. 



2. Notation, assumptions, and some identities 

In this section we provide some notation and background on the renormalized en- 
ergy. We then list some differential identities on the supercurrent and energy density. 



Finally, we discuss some assumptions that we will want our solutions of (1.1) to 
satisfy. 

2.1. Renormalized and approximate energies. Recall 
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and Qp when the center of the baUs is unambiguous from the context. Recall from 
[3] the canonical harmonic map which satisfies the following Hodge system 



div = curl j(u^) = 27r ^ 5, 

with either 



j(n^) • r = ndrO + dr^p 

on dQ or 

j(u^) -1^ = 

on dQ. Furthermore, there exists a G with V-'"G = j{Ui,), where G{x) = G{x;a) is 
defined by the following Poisson equation 

n 

(2.1) AG = 2TTdj in n 

i=i 

with either 

duG = ndrO + dr^ on dQ 

or 

G = on dn. 
The renormalized energy then is defined as 



W{a,d) = lim\ [ - iVu^P — vrnlog - 



We also define the approximate energy as 

We(a,d) = W{a,d) + n(7r|loge| +7) 
where 7 is the constant of [3]. Finally, the excess energy is defined as 
(2.2) D{u;a,d) = Esiu) -Weia,d). 

As u and d are assumed fixed, we usually simplify this to D{a) when the context is 
unambiguous. 

We start with a characterization of the gradient of the renormalized energy. 

Lemma 2.1 (Lin [30], Jerrard-Soner [21]). Let ^ G O"* and d G {±1} then the 
canonical harmonic map u^, = M*(-;^,d) and the renormalized energy W{^,d) satisfy 
(2.3) 



1 2 



^-Y,djdMi3){^iW{i,d))^ 



vr 



where Lp G G^(il) and V^c^ /zas support in a neighborhood of the 's. 

Next we list estimates on the canonical harmonic map and renormalized energy 
and their derivatives. 
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Lemma 2.2. There exists constants C depending only on $7 such that for every 
bounded, open Q C M?, a G il"* and d G {±1}", the renormalized energy W{a,d), 
canonical harmonic map Ui,{-,;a,d) and its potential G{-;a,d) as defined in ( |2.1[ ) 
satisfy 



(2.4) \\jiu*)\\L^(nr{a)) 
for all r < Pa, and 



|VG||icx)(Q^(a)) < 



2n 



(2.5) 



\ViW{a,d)\ < 



Cn 
Pa ' 



Cn 

\ViVjW{a,d)\ < -j- 

Pa 



for every i,j £ {1, . . . ,n}. 
We also have the bound 

(2.6) 



W{a,d) < C{n^ + 



-) 

pV 



Finally, let ^ = {(i, . . . , ^n) and ^' = {^[, . . . , ^'J with ^' £ Jl"*. Let flri^, 
n\ (u'^^^Bri^.j) U Br{^'j)Y then for every d £ {±1}", 



(2.7) 



1 " 



3 = 1 



for all r < mm{p^, p^'} , and additionally, for 1 < p < 2, 

(2.8) ||j(n(0)-iKO)llL.(n) < {7.Y.\^^-C\)'-\2nn: 



2-^ 



Proof. The Neumann boundary condition results are proved in p3j. Corresponding 
results for the Dirichlet boundary condition can be established by using similar argu- 
ments along with estimates of the corresponding renormalized energy found in [47J. 
Further estimates in the Dirichlet case will be discussed in Section IH □ 

2.2. Differential identities. The following identities hold for smooth solutions of 

2 l)=2|Vnp 



(1.1): 



(2.9) 
(2.10) 

(2.11) 



mass identity 
supercurrent identity 



1 



a- A "'^■'2 



|loge| 
1 



\u\ 



u 



\loge\ 



{iu, dtu) = divj(u) 



energy identity dtee{u) = div(nt,Vn) 



\dtu\' 
|loge| 



For fixed e regularity follows from standard parabolic theory. We remark that (2.9 ) 



can be used to show that < < 1, (2.10) will be used to show that j{u) is nearly 



divergence- free in a time- averaged sense, and (2.11) is the primary tool to establish 
the vortex motion law. 
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2.3. Assumptions. We construct initial data that ensures that solutions to (1.1) 



with the following bounds on fundamental quantities. Though values certainly are not 



the optimal choices, they do represent values for which we 


can close in our argument 


(2.12) 


number of vortices 


n < log 


1 

£200 


(2.13) 


minimal intervortex distance 


> |loj 


1 

rgj 100 


(2.14) 


total time scale 


T < 1 lo^ 


5|loge| 1 


(2.15) 


time averaging scale 


6e = \log 


1 


(2.16) 


resolution of vortex location 


Ve < |lo 


1 

ge| 4 


(2.17) 


initial excess energy 


D{am 


2 

< loge 5 



Given these choices, we will establish a few composite estimates in the next section. 

3. Measuring true vortex positions 
We define a series of time intervals on which our function u is well-behaved in 



different senses. In Subsection 3.1 we will define ??(t), which is a natural measure of 
the distance of the vortex positions to the a^-'s. 

(3.1) TO = inf {t>0: p^^t) < 

(3.2) Ti = sup < < t < To such that \\J{u{s)) -y2'^^i^<^^(s)\\w-''-Hn) ^ 

* I i=i |loge|2oo 

and D{a{s)) < 100 for all < s < t| 

(3.3) r2 = sup 1 < t < Ti such that i]{s) < -V^ for all < s < t 

t [ 2 

(3.4) T3 = sup he < t < T2 such that (t?)^^ (s) < for aW 6e< s <t) . 



We note that is introduced in Proposition 5.2 By the end of Section [8] we will 



show that these time intervals overlap, at least until T defined in (1.15). 
The definition of ri implies that 

n 

(3.5) pa^t) >P*> llogepTSo , ||J(u) - J]]7rdi5„^(,)||^_i.i(^) < llogep^ 

i=l 

for all t G [0,ri], when e is sufficiently small. Therefore when this holds there exist 
^{t) = (^i(t), . . . ,^„(t)) e f)"* such that - Oil < ^ for all i, and 



(3.6) 

II ^ 

I loge 



i=l 

n 



ee{u){s) _ ^ 

i=l 
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Ce 



(3.7) 



n 



C 

\\og£\ 



n 



nlog — + W/(a(0)) 
P^ 



Those bounds are found in Proposition 5.2 and Theorem 5.3 



We bound the fohowing composite values based on the time interval definitions. 



Lemma 3.1. Assuming (2.12|)-(|2.17|) then for < t < ti 

(3.8) 
(3.9) 

(3.10) s,<ein, 



W{a{t)) < |loge|ioo , 
Ee{u{t)) < |loge|i+^ , 



(3.11) 



te<|loge| 100. 



Proof. By (2.6), since Pa_(t) > then (3.8) follows from (2.12) and (2.13). Next note 
that 

E{u{t)) = We{a{t)) + D{a{t)) = n (vr |loge| + 7) + W{a{t)) + D{a{t)) 
< Hog el 



from ( |2T2| ), ( |3^ , and t he fact that D{a{t)) < 100 for < t < n. Finally, ( |3l0l ) and 
(STTl) follow from (|2l2|, (plsl), (|3^, and ([3^. □ 



3.1. Definition and properties of r}{t). Since the energy is concentrating at the 
points and the ODE gives us vortex positions Oj, the main objective of our the- 
orem is to estimate and control \^{t) — a{t)\. This is a challenging quantity to work 
with directly, so following Jerrard-Spirn [23], we define a similar quantity that is 
differentiable and has very similar properties. However, in [23] the authors use the 
Jacobian instead of the energy density due to the differing conservation laws for the 
Gross-Pitaevsky equation. In the parabolic problem here we set 



(3.12) 



where 



\\og£\ 



$j(x,t) = ^p{x - aj{t)), ip{x) = xXpAx) 



1 for \x\ < 1 



. The 



and XpAx) = x{fj for a fixed x e C^{^^) satisfying x{x) ^ l q for |x| > 2 

^>j's are supported on B2p_^{aj{t)), so that {supp <I>j(x, t)} are pairwise disjoint when 
Pa{t) ^ P* and in particular for all < t < ri. 
We claim 
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Lemma 3.2. If < t < ti then 

vit)-^T^\m-aiit)\ 

i 

.e{uit)) 



I log el 



i=l 



^{t)-\\J{u{t))-Y,^d,5 

ai{t)llvi/-i.i(Q)| 

i=l 

Proof. First note that in view of, the definition of ri, and 

n 



(3.13) 



4 = 1 



< llogel 200 +U < 



when e is sufficiently small, for all t € [0, ri]. From the definition of $j it follows that 
^j(i) — aj(i) = ^j{^j{t),t) for all such t. Therefore, there exists a unit vector Vj{t) 
such that \(j{t) — aj{t)\ = vj ■ {2dj^j{S,j{t))); hence. 



j 



<r,{t) 



-/('E%,..-^)(E.-, 



ee{u{t)) 



IT 



w- 



|loge| 

(3.14) < r]{t) + Cte for all t e [0, n]. 
A similar argument shows that for such t, 

(3.15) r?(t) <^^|ei(t)-a,(t)|+C7t,. 

Again following the idea in [23] we use the triangle inequality and the norm as 

the length of a minimal connection to get 



Mu{t)) 

|loge| 



i=l 



< 



ee(n(t)) 



i=l 



\loge\ 

< t, + 7TY,Mt)-ai{t)\ 
(3.16) < Cts + 2r]{t) for aU t G [0, n]. 



1=1 
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In the same way one finds that 

(3.17) vit) < Cte + ll^n^ - E^^^'^a,(t)llty-i,i(o) for ah t G [0,ri]. 

In particular this implies that when e is sufficiently small, 

(3.18) rj{t) < 2Ve for ah t £ [0,ri]. 

Finally we can follow a similar argument to show the last estimate. □ 

4. Initialization of Gronwall argument 
Use the notation r]j{t) := J ~^^^^^j{x,t) £ as introduced above, then we have 



Lemma 4.1. Let u be a solution to ( 1.1 ). Then for < t < ti and j = 1, . . . ,n 

(4- 1 ) rj3= Tj,i + Tj- 2 + T,- 3 + T,- 4 + + Tj, e + Tjj 

where 

= y^iij - aj) ■ J {VjWiO - VjW{a)) 

7;-,3 = J d^^^^j {d,,\u\,d,M)-dl,,^MH) 



3^^) ■( N 



0,5= (j(^.)). 



Proof. Differentiating r]j, we obtain 

d _ f ee{u) d I' d ee{u) 

dt''' ~ J ]logi\ dt' + J 'dt\logi\ 

Since ^$j(x,t) = ^f{x — aj) = {—dj) ■ Vf{x — aj), we can use the ODE and the 
fact that ^j{^i{t)) = for i / j to write 

^eiu) d ^ f egiu) , . , „ , , 



\loge\ dt J |loge 

= —djVjW{a) ■ Vip{(,j — Oj 



+ / (^-E-'^e.) (-^V,ma)) •V^(x-a,) 
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Next from the evolution identity for the energy, the representation Vn = ^V|u| + ^"^j^^ 
and -B"^ = \A- + 2{A-B)-B we find 



= j Cx.^j- 



{dtu,dx^u) 



'^^ \\oge\ J^' 



|loge| 



{d^i,u,dx^u) - ^ |Vn|^ 



o2 (l-kP)' 



|loge| 



+ 



1^1 V l^^l J k 2 



1 " 

fc=i 



□ 



Next we obtain an estimate of ?) by separately considering contributions from the 
different terms isolated in Lemma |4.1[ We will prove 

Lemma 4.2. For t € [0, ri] 

I'iWI < — Me sup ll{s)+Be \+—lAe SUp ri{s) + Be- 

P* \ se[o,t] / P*y se[Q,i\ 



Proof. The condition e < £9 is needed to guarantee the validity of estimates (3.14), 
(SlTl), delol) from Section [tl 



Note from Lemma |4.1| and the definition (3.12) of 7/ that 

n 

(4.2) ?) = Ti + . . . Ty, where T^fc = ^ • 

We estimate these terms in turn. 
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First, note that V(j){^j — aj) = — aj for < t < ri, by the definition of (p and 
( 3l3l ). Thus, in view of ^J^ , 



a . 



And arguing as in the proof of (2.7) we see that 

n 



|V,-M^(0 - VjW{a)\ < V Mt) - afc(t)|(sup sup |VfcVjVF(y)|) 

fc |y-aW|<|€W-a(t)| 



< (r?(t) + Ct,)C 



n 

72t 



using (3.14) again, as well as bounds on V W from (2.5). Thus 
(4.3) 
Next, 



\Ti\<C^{ri{t) + Ct, 



J{u) -Y.^d,5^^ ] I Y.^V,W{a) ■ V(cl>, • 



< 



vj;V,H^(a)-V(ci>, 



L°° 



Since the <I>j's have disjoint support 



Vj;V,W^(a).V($,.-^) 



< sup \VjW{a)\ llV^^jlloo < C 



n 



We conclude from (3.6) and the above that 
(4.4) 



n 



\T2\ < CSe^. 

pi 

Continuing, we use the fact that V^$j vanishes in Bp^{aj), together with (6.19), to 
find that 



\T3\< 



(4.5) 



I 'b I 



\V\u\ 



flp^ (a) 



< 



1 



As sup r]{s) + Be 
se[o,t] 



Exactly the same considerations show that 

1 



(4.6) 



IT4I < 



Ae sup 7]{s) + Be 
se[0,t] 
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Next, 



\n\ < 



3L 



IIjK)II 



L2(u^.suppV2$j)- 



Using (2.4), one can easily check that ||j(if*)||L2(UjsuppV2<i>j) < ^{Cnp^)2 , and hence 
we conclude that 



(4.7) 



\n\ < 



3 

n2 



As sup rj{s) + Be 

se[o,t] 

3 

n2 



Exactly the same argument shows that \Tq\ + jTyl < ^JAe sup^gjo,*] '^(•5) + 



□ 



The result of the lemma is not good enough to get any very strong result from 
Gronwall's inequality, but it still implies a useful Lipschitz bound. Since we will be 
using a time averaging operator, we define the time average of a function h as 



{h)s^ it) 



1 



h{s) 



Jt-S, 

for any t > 6s. This will be used to handle estimates that make use of the differential 
identity (2.10) to show that the current j{u) is effectively divergence- free on average, 
and such estimates. 

Corollary 4.3. We have for t G [0, T2] 

3 

712 , 7 1 

(4.8) < —yjAsVs + Bs < llogepso |iog|loge| |2. 
IfO<t-6s<s<t<T2 then 

3 

n2 I 27 1 

(4.9) \7]{s) - {r]{t))s^ I < 5s — VAv7+Bs < \loge\-^> \ log |loge| 1 2 



Proof. Both estimates come directly from (2.12), (2.13), (2.15), Lemma 4.2, and The- 
orem |6]T1 



□ 



5. Quantitative rate of second-order gamma convergence 



In order to quantify the closeness of our heat equation solution to the point vortex 
model, we need some results that quantify the closeness of u to the canonical harmonic 
map Ui,, based on the excess energy 



D{a) 



es{u) - Wsia,d). 



We also define the excess energy in the ball Ba{aj) to be 



D 



es{u)dx — (Trlog — h 7 



a 



where 7 is the constant from [3]. The following result lets us compute the rate of 
convergence of u to the canonical harmonic function, based on the excess energy. 



VORTEX LIQUIDS 



21 



Proposition 5.1 (Jerrard-Spirn [23]). Let Q be a bounded, open simply connected 
subset ofM? with boundary. Then there exists constants C, Ki depending only on 
f2 such that for any u G if ^(17; C), if there exist n > 0, finite, with a = (ai, . . . , an) G 
17"* and d G {±}" such that 



(5.1) 



J{u) — irdjdo 



for some Se G [e^/log(jhJe), pa/ Ki], 



and if Ase<a* = J^ (s, + eE,{u)) < th 



en 



(5.2) / ee{\u\) + \ 



\u\ 



f 



<D{a) + Ci —{se + eEe{u)) 

V Po^ 



for canonical harmonic map u^. Finally, 

(5.3) IU(«) - j(^.)llL4/3(f7) < C/DR + error 



1 3 

error < C£2E^(u)^ 



+ Cn {se + eEeiu))i 



n \ 4 1 / 

+ 1 + 



i^.(n) 



We also have the fohowing the result that lets us localize the Jacobian about delta 
functions. 

Proposition 5.2 (Jerrard-Spirn [23]). Let be a bounded, open, simply connected 
subset o/M^ with boundary. Then there exists constants C and C^, depending on 
diam(r2), with the following property: 

For any u G H^{Q,; C), if there exist n > 0, a = (ai, . . . , an) G f^"* and d G {±1}'^ 
such that 



(5.4) 



J{u) — Tldjdo 



< 



Po 



and if in addition E^i^u) > 1 and 
(5.5) 



^ 1 

— n + — \/^ 

Pa 



< 



then there exist (^i, . . . ,^d) G f^"* smc/i i/iat — ai\ < 2^1? "'^'^ 



(5.6) 



i=l 

< C e 



n 



(C + Z)(a))2e^^(") + (C + £>(«))— + Ee 

Pa 



where D{a) = D{u, a, d) = E^(u) — We^a, d). 
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We now state a result that clarifies the convergence of to a set of delta 

|log£| 

functions. 



Theorem 5.3. Let u satisfy the same hypotheses as Proposition 5.2. then the 
found in Proposition 5.2 satisfy 

P _ / 7/ 1 

(5.7) 



|loge| 



n 



< 



C 



llogel 



n 



n{D{a) + C + log — ) + W{a, d) 

Pa 



Furthermore, for all t G [0, ri] then we can estimate 
ee{u) 



(5.8) 



Hog el 



n 



< 



n 



|loge|/)^ 



< loge 100 



We make precise (and quantitative) an argument found in [8]. The first step is a 
moment estimate on the Ginzburg-Landau energy about the vortex core. 



Lemma 5.4. // \\J{u) — vrtJi 



kr ee(n) -vrln^ 



Kq then there 



exists ^ G -6^/2(0) (^'^d a constant C , independent of Kq and u, such that 



(5.9) 



/ \x-i\ee{u)<rC{K^ + l) 

J Br 



Proof. By Theorem 1.2' of ^22j then there exists ^ G i?r/2(0) such that for any r < 
r — 1^1 and e < o" < r, 



(5.10) 



Ce ("U) > TT In Kq — C 



Prom (5.10) we see that 
(5.11) 
so that 

ee{u) + 



I Br 



(5.12) 



(li) > vr In Kq — C 



ee{u) < -kIh- + Kq - -kItl- + Kq + C 

Br(0)\Br{£.) e £ 

< + CKq + C 

T 
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Now we look at the energy in the annular set i?2-j\-B2-(j+i) (0 ^r\-Se(^). In 
particular 

/ ee{u) = / ee{u) - / ee{u) - / es{u) 

iB2_,A-B2-(i+i)«) J Br JBr\B^.j{0 ^B^-t.+D^) 

< vrin - + - vr In + C{Ko + 1) - vr In + c{Ko + 1) 

= 7rln2 + 2C(Ko + l) = C{Ko + l) 
Next we prove the claim. If we let 2~^^ = | and 2~^= = e then 

/ \x-(.\ee{u)= \x-^\ee{u)+ \x-^\ee{u) 

J Br JB, JBr\Br/2{(,) 

+ \x-C\e,{u) 

r/2 ^-^ 

< rCiKo + 1) 

since E^M, 2"^ < r. 



□ 



In order to establish the proof of the theorem, we use the following energy lower 
bound: 

Lemma 5.5 (Jerrard-Spirn, [22]). There exists an absolute constant C > such that 
if u ^ H^{Ba-) satisfies 

||J(n) ±7r5o||vy-i,i(i3„) < ^ 

then 

(5.13) 0<Db^+ c^^Jh^+ ^ \\Jiu) - vr<5o||^_i.i(5^) . 

We now present the 



Proof of Theorem 5.3. From the proof of Proposition 5.2 in [23], the is chosen so 
that jj- <\. Then choosing a = ^4^, we find that 

(5.14) Ase = -^ = -<a< 



2Ci,n^ n uKi 
In particular this implies 

for a defined above. 
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1. Given the choice of a we claim that the foUowing bounds hold 

C 



(5.16) 
(5.17) 



n 



<DB„(^^)<D{a) + C 



[ ee{u) <D{a) + Wia,d) + C. 



In order to prove (5.16) we note that (5.15) and Lemma 5.5 imply < D^j + 

Cf + ^ II J(u) + 7Tdj6o,^ ^ + n ' ^^lich folloWS if 

a > r?e. 

To prove the upper bound we use the following inequality that can be found in the 
proof of Theorem 3 in 



(5.18) 



„ n 

D{a)= / [e,(n)-e,(^.)] + 5;i)s^(„^) + 0(( — )2) 



> 



eei\u\) + 



1 



3{u) 



if cr < ^. Thus E"=i^Ba(a«) ^ ^(«) + and hence Ds^^a,) < D {a) + C since 
DB„{a ) > for each j G {1, . . . This finishes the proof of ( 5.16| ). 

For'dlTTi) we again write Z)(a) = 6^(7^) - e,(n^) + E 2^5,(^0 + 0((^)^), so 

/ e,{u) = D{a) - DB^ia,) + / e,(n,) + 0((- 



= D{a) - + W{a, d) + 0((— )2) 

Pa 

< D{a) + W{a,d) + C. 



2. From Lemma |5.4l (5.16), and (5.15) there exists a in B^i2{oij) for each 
j = {1, . . . , n} such that 



(5.19) 



/ \x-ij\ 



ee{u) < Ca {D{a) + C) 



Next we choose an arbitrary cj) G M^Q'°°(ri) with || V(/)||ioo(Q) < 1. Then 



< 



|loge| 

, ee{u) 



\loge\ 



eeju) 
|loge| 



A + B. 
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We first handle A. Choose any xq G di), then from (5.17) 



(t){x) 



\loge\ 



{4>{x) - 4>{xo)) 



|loge| 



< / \X — Xq\ 



eeju) 
|loge| 



(5.20) 



^dian^ 
|loge| 



Next we estimate B. Without loss of generality assume aj = and = ^ and 
again choosing xq G dQ, then 



eeju) 
|loge| 



< 



{cp{x)-m) 



|loge| 



+ \m\ 



eeju) 
|loge| 



< 



\x-^\ 



ee{u) 



|loge| 



+ \^-xo\ 



Ba 



\loge\ 



From (5.19) we have 



|loge| 



Whereas, 



/^^ ee(M) - vr |log e| < vr log ^ + Db^ + 7 imphes 



Bo < 



diam(r2) 
Hog el 



Since a = rf" a < 1 then 



Bi + B2< 



C 



|loge| 



log - + D{a) + C 
a 



log — + D(a) + C 



Po 



Combining with the bound on A yields 



ee[u) 
Hog el 



< 



C 



|loge| 



n 



n{D{a) + C + log — ) + W{a, d) 

Pa 



and since n < — then (15.71) follows. 

— pa vi^^jjj; 



□ 



6. QUANTITIVE BOUNDS ON THE KINETIC AND EXCESS ENERCY 

We now present a kinetic energy bound for fixed e. Similar bounds with errors of 
the form Oe(l) can be found in [181 131} 

Our method of proof is inspired by the choice of test function found in the proof 
of Theorem 3 in [25] . 
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Theorem 6.1. Fix 6 > and t < ti of (3.2), then for any time interval ti < t < 
t2<T with \t2-ti\ > 1, 

H2 



TT 



(6.1) 



where 



(6.2) 



Id|2< 



ti 



/ X[x-a{t))- 



+ CAe sup T\Cjit)-aj{t)\+CBe 



Be 



< lloeehs lloellogel 



and where C depends only on 0,. 
Furthermore, if D{a{0)) < Be, then 



T ^ |loge| 50 |log|loge| 



(6.3) 

and 



D{a{t))<Ae sup r]{s)+Be, 

s6[0,t] 



(ie{\u\) + 



< Ae sup r/(s) + Be 

s6[0,t] 



(^) - J'(^*)IIl4/3(Q) ^ \Ae sup 7/(s) + 

V ^e[o,t] 



A similar theorem was proved in j28] for a single vortex that stays an 0(1) distance 
from the boundary for an 0(1) time. Here we prove a much more exphcit estimate. 
The major tool to estabhshing a finite-e bound on the kinetic energy is the following 
optimal result on the equipartitioning of Ginzburg-Landau energy, which improves 
related results in [3T] and 



Proposition 6.2 (Kurzke-Spirn 

{u) < IT log J + Kq then 



I. Suppose \\J{u) — ttSqI 



(6.4) 



\dxiu\ 



I ((9a;iM, 9^2^)1 \dx2u\ 



flogf 



flogf 



< Wi^ilo, 



where Ki = C{C + Ko)e'^o/n g^^d C is a universal constant. 

We apply this equipartitioning result to the evolution identity for the energy and 
deduce a rate of convergence for the kinetic energy. 



Proof of Theorem 6.1. 1. We generally assume t < ti. Per (3.2), we have for every 
s € [0, t] the estimate 



(6.5) 



< a 
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Conservation of energy implies 



(6.6) 



Jn 



\dtu\' 
\\oge\ 



D{a{0)) + We{ao) - We{a{t)) - D{a{t)) 



From Q we see that D{a{t)) > -C (§e + eEe{u{t))) > -C. We also have 



from g that W{ao) - W{a{t)) < \W{aQ)\ + \W{a{t))\ < C(n^ + ^). Finally 
D{a{Q)) < Cn^ < C. Since n < yO~\ we find that 



(6.7) 




2 2 

< 



\dtu\ ^ n 



Jo. 



7r|loge| ^ pI' 

and from (5.17) and the assumptions on Z)(a(0)) and W{a{s),d), we have that 

,2 



(6.8) 



ee(n) 



< 



n 



Pi 



for every a with a > C^. For a < p^,, B2a{aj{s)) n B2a{ak{s)) = for all s G [0,t] 
unless j = k, and i?2o-(aj(s)) H c^fi = for all j. 

2. We now make the following claim. Let x ^ C^{M?) be a function such that 
X > 0, X = 1 on 5^(0), X = in \ 52r(0) and \D''x\ < Cr^^ for k = 1,2. Then 
there exists a constant such that for any hj G C^([0, T]; M^), j = 1, . . . , n: 

r j Y.x{x- ajit)) bjit) ■ (dtu, Vn) + vr r bjit) ■ dj{t) 



< 



sup \\b 



T r oo r oc 



a 



(6.9) 



+ te 



a 



- sup ||aj||^ooj;^oo sup ||6j||^ooj;^oo + SUp 



r CX3 r oo 



+ Tsup 

j 



1 

|loge| p2 



ll^.l 



T oc r oo 



+ 



,1/2 



sup 



r oo r oo 



pl(T\\oge\ ' 3 

For any test function (j) G C^([0,T] x 17) with compact support in $7 and any 



< ti <t2<Twe have 



(6.10) 



dt 



ti Jn 



I logs 



ti 



|loge| 



1 



|loge| 

^2 



ti Jn 



|loge| |loge| 



as is easily seen by multiplying (2.11) by j]^^^ and integrating by parts. 
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We now follow [25| and set 



<Pit,x) = ^(pjit,x) = ^x{x- aj{t)) bj{t) ■ {x-aj{t)). 
Then we calculate, dropping the t-dependence of a and b, 

n 

V(p{t, x) = ■ (x - aj))Vxix - aj) + bj xix - aj) 

n 

dt(t){t, x) = ^ -hj ■ Vx{x - aj) bj 



+ x{x — CLj) ( bj ■ {x — aj) — bj ■ cij 



Vdt4>(t, x) = — dj • V^x{x — aj) bj ■ {x — aj) — hj ■ Vx{x — aj)bj 



+ Vx{x — aj) ybj ■ {x — aj) — bj ■ hjj + x{x — aj)bj. 

In particular we have (with constants depending on a and b) 

1 



IIWIIl-l- ^ -sup||6j||^oo^c> 



1 



1 



r oo r oo 



Now we analyze the terms in (6.10) one by one. We have by (|3.6| 




dt 



< 



h Jn 

Ht2 



eeiu) 



+ 



1) ■) 



IIWIIls?l- sup 



llogel 
eeHt)) 



|loge 

ee{u{t2,-)) _ 
|loge| 
ee{u{ti, ■)) 



te[ti,t2] 



|loge| 



w- 



a 



sup \\bj 



On the other hand 



Y,bj{t2) ■ iCj{t2) - aj{t2)) - bjih) ■ (e,(ti) - ajih)) 
i=i 



< 



sup\\bj\\^^^oo sup V|0(i)-ai(i)l- 



Therefore, 
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(6.11) 



rt2 




dt 



h Jn 



|loge| 



< 



sup \\bj\ 



r cxD r oo 



For the second term on the left-hand side of (6.10), 



and 



Thus 



Jti Jn I log e I Jti Jn 



<T\\Vdt(l)\\Lc^L?? sup 



te[ti,t2] 



<trT 



|loge| 
sup \\bj 



w- 



1 „• „ 1,1, „ 1 

"2 sup \\aj\\L^L?f sup \\bj\\L^L?f + " sup 



t2 r " rt2 

dtcj) ^ vr5g^. + 71" ^ / &i • a 



ti Jn 



t2 



^ ^sup 
j 



r^roo sup V|e,(t)-a,-(t)|. 



(6.12) 



Jti Jn |loge| ^Jt^ 



T 



sup ||aj||^cx)^oo sup ||6j||^cx)^oo + sup 
a j T n j T n j 



+ Tsup 

j 



Note that the previous equahty contains the term we want to appear on the left-hand 



side of (6.9) 



For the first term on the right-hand side of (6.10) we use (6.7) and get 

\dtu\'' 



(6.13) 



1 



r.i2 




< 



1 



|loge| Jt^ Jn |loge| ~ |loge| 



I - 
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Finally, for the second term on the right-hand side of (6.10) we have 



(6.14) 



1 



|loge| 



Hi Jn 

n ^ 



V<j) ■ {dtu, Vu) 

t2 




Vx(a: — aj) ■ {dtu, Vu)bj ■ {x — aj) 



^t2 



^ |loge| Jt^ Jn 



X(x - aj)bj ■ {dtu,Vu). 



The second term on the right-hand side of (6.14) is what we want for (6.9). We 
estimate the other term. Since 



1 

p[ |loge| 
1 



^t2 



*1 B2r{a^(t)) 



Vx(x — Oj) • [dtu, Vu) bj ■ {x — aj) 



< 



1 1/2 



sup 



|loge| ' J 

H2 



j 



E 



^ \Jtl J B2r{aj{t))\Br(aj{t)) 



|loge| 



1/2 



(■t2 



1/2 



\Vu\ 



< 



< 



< 



c7|loge 
1 



1/2 



1 1/2 



sup 



sup 



rt2 




ti Jn 



\dtu\^ 
|loge| 



Itl J B2r{aj(t))\Br{aj(t)) 
'■t2 



+ 




ii Jna{aj{t)) 



|Vn| 



o-|loge| ■ 3 
n^\t2-ti\ 



n n 

^ H — 5" ^2 — ii 

Pi Pi 



sup 



r oo r oo 



3. We now study the momentum term on the left hand side of (6.9). From Cauchy- 
Schwarz 



|loge| Jt^ Jnfri 



X — Uj) bj ■ {dtu, Vu) 



< 



rt2 



|loge| 




h Jn 



Xix - a)\dtu\' 



1/2 



rt2 



\loge\ 




ti Jn 



x{x - a)(6 (g) 6) : {Vu Vn) 



1/2 



where (6 (g) 6)jj = for b £ M? and (Vu (g) Vu) - = {diU, dju) for u £ C For any 
6j G and x above, we claim that 



(6.15) 



/ / Vx(a;-ai)(6i^?^i) : (Vn^Vu)- / V^|& 



I log 



^ "^,,, sup||6j||^ 
|loge| ' 3 ^ " 
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2\dyU\ 



+ E 

^ /l + /2 



First we analyze Ii. From (6.4) and choosing cr = S> e, 

^1 ^ e / - (^-^' ^^^) 

|10ge| f-f JB,„fa,-) 



i32,.(a,)\i?<.fe-) 



X(x - aj) 



|loge| |loge| 



\b^\'\dM' ^ \b'j\^\dyu\ 



\loge\ |loge| 



< — 
~ |loge| J 



sup II 6 



■JllL5?L^ A/ log 

r-i2 



+ sup||&j|lLopic»y^^ yn^aW) |loge| 



< 



sup ||6j| 



-t^B.fe) log? 



n n 

+ 



+ 



nloe 



|loge|^/^ P*|loge| |loge 



Next we look at l2- Again from (6.4) and since a ^ e, 



C7\-l/2 



III, l|2 

n7^sup||6,-||^^ 

|loge| ' J ^ " 



Comparing the terms from Ii and I2 yields (6.15). 

4. We now use (6.9) and (6.15) to derive the expected kinetic energy estimate. 
Setting bj = aj, noting that aj = -^Va^W and bj = dj = ^VM^V^iy, then ( [275 ) 
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II . II < ^ 

sup WO'jULZ'L^ ~ „ 
j ^ " 

sup \\aj\\i^<^j^a^ ^ —r^. 



So using (6.9) and (6.15) we have 
t2 / rt2 



TT 



, • |2 

a 



'1 



< 



ti Jn 



X{x - a) 



where 



Fi = C 



C 



p6|loge|l/2 



+ C 



\\oge\ 



sup 



2\ 1/2 



7r^*'|6|2 + F2j 



1/2 



* se[ti,t2] j 



,(s) - aj(s 



p2|log£|l/2- 

Square the previous inequahty, obtaining by division 



< 



|loge| 



t2 



ti 



X(x - a)\dtu{' 



Setting K = |ap and assuming |dp > C > 0, we have using 

K + F, -K + F,-2{F, + F,)+ ^ ^ 
that actually (dropping the final term in the previous equality) 



1 



|loge| 



t2 



x{x — a)\dtu\^ > vr 



CFi - CF2, 



and so (6.1) follows, since Fi> F2. 

5. We next will show that u{t) is well-approximated in certain ways by the canonical 
harmonic map Ui,{t) := ■u^(-; ^(t), d) for t < ti. To do this, we need to estimate the 
surplus energy D(^(t)) with respect to the points ^{t) found in Step 1 above. 

Fix t £ [0, Ti] and assuming D{a{0)) < Be then 

Dim) = D{a{t)) + W{a{t), d) - W{m.d) 

< CAe sup r/(s) + CBe 
se[o,t] 

+ 1 (sup sup \Vy^W{y)\). 



where 



J \y~a(t)\<m-<t)\ 
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From if y G $7" is such that \y — a{t)\ < \£^{t) — a{t)\ , then py > ^Pa(t) ^or all sufficiently 
small e, and so we get 



(6.18) 



Di^t)) < CAe sup 77(s) + CBe + {Cte + r]{t)) — 
se[o,t] P* 

< CAe sup ri{s) + CBe. 



Thus the first conclusion 



/ 

Jn 



eeiHt)\)+- 



\u{t)\ 



n 



(6.19) 



< CAe sup r/(s) + CBe + C[—{se + eEe{u)) 

se[o,t] VP* 

< CAe sup r/(s) + CBe 

se[o,t] 



for all i G [0, ri], where (since > CeEe{u), see the definition (3.7)) (Note that a 
condition a* < pa appearing as a hypothesis for conclusion as a result of the definitions 
of Ti, Se etc.) 

From the other conclusion we deduce that 



(6.20) ||J-(^^)W -iK)WllL4/3(Q) <C\Ae sup r,{s) +C^Be. 



□ 



7. Quantitative bounds on the supercurrent 



In this section we prove estimates of Ti — Tj after averaging in time. As in |23| 
we cannot rely on the Gamma convergence estimates to bound these terms since the 
square root in rj pops out. Following p3] we establish direct bounds on j{u)—j{Ui,) via 
a Hodge decomposition and use a time-averaging. Unfortunately, since the argument 
in [23] relies on very different assumptions, we need to reprove several bounds with 
the worse behavior in the parabolic problem. Our primary result in this section is the 
following bound on all Tj's: 

Proposition 7.1. For all t G [(^e, ti] and j = {1, . . . , 7} 



(7.1) 



{T,),^{t)\<Ae sup (r?),Js) + |loge| 
se[<5,i] 



11 

' 40 



where Ae = which satisfy (2.12), (2.15), (2.14). 



Terms Ti through T4 can be directly estimated from the non-averaged estimates in 
Section [4j However, we are unable to bound the rest of the terms without resorting 
to explicit time-averaging estimates of the differential identity (2.10). The following 
lemma addresses the much harder terms. 
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Lemma 7.2. For all t G ri] and j = {5, 6, 7} 



(7.2) 



{Tj)sAt)\<A sup {v)s As) + \loge\ 

s&[S,t] 



11 

^ 40 



w/iere A = i(;/iic/i sate/y p32| ), I^J§, ( |234| ). 



Proof. For simplicity we write 

where 
(7.3) 



j(^^*) 



= 1,2, 



Ji2 = — JI21 = 1, Jii = J22 = 0, and jm denotes the m component of ji^u^,), m = 1,2. 
Here n*(x,t) = M*(x;^(t)) as usual. 

1. From the definitions and (2.4) one finds \0 < C-^, and IsuppO < Cnpl- It 
follows that 



(7.4) 



\Li{n) < n^'P* 



for 1 < q < 00. 

The following proof is quite similar to the proof found in Proposition 1 in [23]; 
however, we include it since the bounds are different, due to a different conservation 
law. We perform a Hodge decomposition 

(7.5) j{u)-j{u.) = Vh+V^f2 
with boundary conditions either 

(7.6) /i = and 8^/2 = on dn, 



or 



(7.7) 



d^fi = and /2 = on dQ, 



depending whether we are dealing with Dirichlet or Neumann boundary conditions. 
And so we examine 



A/i = divj(^i) 
-A/2 = 2 \j{u)-Y,T^djSi:^ 



with (7.6) or (7.7). 
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in), 



3{u) 



5. 



il-\u\)) +( / C-curl/2 ) +( / C-V/i 



5. 



/C-^^(1H.|)) +(/C-curl/.) +(/(C),.(V/.),^ 



+ (/ (C-(C)4)-(v/i-(v/i), 

= Ai+A2 + A3 + A4. 
The first term is estimated by the Cauchy-Schwarz inequality, 



(7.8) l^il < IICIU- 
2. Next we claim that 



n 



L2 



Pi 



(7.9) 



\A2\ <CSe 



715 p^^ {Ee{u) + Wk)'- 



< 



\loge\ 



From the Hodge decomposition and standard elliptic estimates [H] we have 



II curl /2 II LP (n) < ||/2||H^i,p(n) 



for 1 < p < 2, with a constant depending on p. Taking ^ = 1 - ^ in (0| for p £ [1, 2) 
to be selected, we conclude that 



|^2| < ||Cl|L^I|curl/2|Up <Cn'"^ s,!-^ {E,{u) + mrf-l 



Choosing p = |, we arrive at (7.9). 



3. Next, we estimate ^3, and here we fundamentally use the time-averaging to 
control V/i. 



|A(/i) 



|div {j{u) - j{u^))^^ 
1 /■* {iu, dtu) 



Il2 



|(divj(n))c 



Il2 



< 



Jt-5, |loge 

1 



L2 

* l^tnl' ^ 1 
4|loge| Jt^s, |loge| ^ 5e\^oge\ pi' 



By standard elliptic estimates 



(/l)5J|£72 ^ ||A (/l)<5j|r2 ^ — 



n 



51 p^ \\oge 
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Combining with (|2J2|), (|233|), (|2J5]) yields 

<c^(IICIIl4/3)5J|(v/i),^ 



< 



11 

n i 



^13 



1 ^ 

I 2 



_ 277 

< |loge| 800 



6ipi I log e I 
hence, 

(7.10) \As\<\loge\-'^ . 

4. Finally, we consider the challenging term A4, and we again following the strategy 
of [23] . The idea is to take advantage of the fact that (5e is small to show that C is 
close to (C)^! and similarly V/i and (V/i)^^. First we have 

\A,\< sup \\as)-{OsJ\LHn) sup ||V(/i(s) - (/i),^ 



(7.11) 



< sup ||C(5)-C(s')llL*(n) sup ||V(/i(s)-/i(/))|| 4 



In estimating the quantities in (7.11), we will use that for s,s' £ [t — 6£,t] with 
t e [Ss,Ti], 



(7.12) 



n 



\aj{s) -aj{s')\ < ^6e, 



which follows from (1.11). Note that (7.12) and (2.5) imply that \dj\ < From 

( |7.12[ ) and ( |3.14D , ( |3.18D it fohows that for s, s' as above, 

"2 2 
(7.13) \2\^.(s) - Cj{s')\ <C-6e + 7j{s) + r]{s') + Cte<-S,+Ve. 

4a. We estimate ||V(/i(s) — /i(s'))ll ^• Assume that s,s' £ [t — S£,t] for t £ 



[(5£,Ti]. By elliptic regularity, ( |7.5[ ), and either (7.6) or (7.7) we find that, 

C7||A(/i(.)-/i(.'))||^, 
\\V-[j{u){s)-j{u){s') 



l|V(/i(.) - /i(«'))ll,|(,) < ^^l|A(/,(.) - /i(^'))ll^-i,|(,) 



(7.14) 



<\\j{u){s)-j(u){s'))\\^.^^^. 



Using the triangle inequality and (6.20), it follows that 

\\3{u){s)-3{u){s'))\\.^^^^ < ClAe sup r/(r) + Cv^ 



'L3(n) 



re [(5,s] 



+ ||j(n.)(s)-j(^.)(.'))ll^4^^^ 



< |loge| 200 |log|loge| I 
+ \\j{u*){s) -j{u^,){s' 



4 
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The last term on the right-hand side can be estimated by combining (2.8) and (7.13), 
and we get 



(n.)(s) - j(n.)(s'))ll^4^^^ < nt ( 5,- 



n 



3 

712 



< —6i +n-2Vi < |loge| 



The rest of the terms on the right-hand side of (7.14) are smaller using the bounds 
on n, Therefore, we find that 



(7.15) 



||V(/i(s) - /i(5'))ll^|(^) < lloger^ I log|loge| 



4b. We estimate ||C(s) — C('S')llL4(r2)- Assume that < t — (5^ < s, s' < t < ri. In 
order to find a time-Lipschitz bound on we have from the definition (7.3) that 



j 



= Z1+Z2. 

First consider Zi. From the definitions. 



^. ($,(«) -$,(.')) 



L°°{n) 



< \\dxix,r.[fix - aj{s)) - ip{x - aj{s'))]\ 



L°°(r2) 



< CWdxix^x^fhoo \aj{s) - aj{s')\ < -^6e 

P-k 



using (7.12). 



As in |23 
(7.16) 



we claim that 

supp V2cD^.(s) U supp V2ci>,(s') C i?3p.(^j(s)) \ {i,{s)) 



for all e sufficiently small. This follows from (7.12), (7.13), and (3.14). In partic- 
ular the distances separating ai(s), aj(s'), ^i(s), ^i(s') are significantly smaller than 

|logep2oo < p^. 

The support condition (7.16) implies that |j(ti^)(^(s))| < ^ on the support of Zi. 
Since the support of Zi has measure bounded by Cnpl, we conclude that 

(7.17) < — [CnplYS, < —6e < |loge| 



163 
'800 



p: 



pi 
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Finally we consider Z^- Since 



< , and using that supp Z2 has 



measure at most Cnpl, we use Holder's inequality to estimate 



2|Il4 



< 



c 



L°°(UjSuppV2<I>j( 



It then follows that supp Uj V^<I>j(s') C f^p^/2(?(s)) n f^p^/2(C(s'))- We therefore use 
([2^ to find that 



\j{u^){s) - j{u^){s)\ 



L°°(UjSuppV2$j(s')) 



Consequently, (7.13) and (4.9) imply that 

11^211^4 < C{^)H-6, + r,{s)+ri{s') + Ct, 



5 



< sup (r/(t))5 + ^ 

,2 8<s<t „2 



(7.18) 



< 



pi 



n 
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_ sup (r7(t))r + |loge| iso | log |loge| 1 2 . 

s<s<t 



Combining (7.17|) and (|7.18|) yields 
(7.19) 



C(s) -C(^')IL4(f,) < C-^sup(??(t)) +C7|log6|-8oo . 



5. Finally we combine the above with (7.11) and (7.19) to deduce that 



\Aa\ < 



~ Pi 



sup {ri{t)) + |loge| 10 



which completes the proof of Lemma 7.2 



□ 



Finally we can finish the proof of Proposition |7.1| which entails showing that time- 
averaging does not damage estimates from Section |4j 



Proof of Proposition 7.1. Since (Ts)^^ through (TV)^^ are covered in Lemma 
concentrate on turning the Ti — T4 estimates into a time-averaged estimate. 
Since tj is continuous, we have for some c S 



7.2 



we 



—As sup r]{s) = —Aeri{c) 

P* se[SeA P* 



<-Ae 
P* 

^ Pt 



(c) + C |loge| 80 |iog|loge| |2 



sup (77)5 (s) -I- C |loge| 80 I log |loge| 1 2 



2rp _27 1 _11 

and since |loge| 80 | log |loge| 1 2 < |loge| 40 ^ the result follows. 



□ 



VORTEX LIQUIDS 



39 



8. Completion of Gronwall argument 



We now complete the proof of our Theorem The proof entails an argument 
that extends T2 up to time T which ensures that the required estimates hold. 

Proof of Theorem \l.l\ Let T = minly^log |loge|^, tq} denote the longest possible 
time interval of existence. The main point of the proof will be to show that we can 
extend the time up to T by a combination of continuity arguments for the Jacobian 
and a Gronwall estimate on {rj)^^. 

1. We first claim that (1.1) is a continuous operator from 1 1— t- H^, i.e. 
(8.1) \\Vuit) - Vu{s)\\l2 < CeOj{\t - s\) 

for < s < t < T, where can depend on e and T. It is standard theory 
(see for instance [H] Section 5.9, Theorem 4) that if u G L'^{0,T; H'^{il.;C)) and 
dtu £ L2(0,T;L2(Q;C)) then u G C°{[0,T]; H\n;C)), which implies These 
conditions are true for solutions of (1.1) since by the gradient flow property, 



r [ \dtu\^ < CE,{u{0)) < Ce 
Jo Jn 



and from (2.9) one finds ||u||loo < 1 due to the maximum principle; hence, 
rT 



Jn 



lAuP < C 



T 



Jn 



\dtu\' + C 



T 







1 , , 

n £ 



:(1 



|2\2 



< CEeium + CT\\u{t)\\l^^^^E,{uo) < Ce 



where Ce can depend on T and e. 

2. Since J{u) = detVtt = —V'^ui ■ Vii2, where u = ui + iu2- Then for any 

< s < t < r 



\\j{um-j{u){s)\\ 



sup 



(«)- 



<i 



ct>{J{u){t)-J{u){s)) 



< C \\Vu{t) - Vn(s)||^2 \\Vu{t) + Vu(s)||^2 

< CeUj{\t-s\) 

with modulus of continuity uj{s) from above and Ce, a constant depending on e and 
T. 

3. Suppose that r2 = ti < T, then we claim that we can extend ri 1— ri + /i for 
some n > 0. Suppose not, then at i = ri we have 



(8.2) 
and 

(8.3) 



Diain)) < 100 



J(u)(ti) -^dj5a^ 



(n) 



< epilog e| 100, 



w- 



where is defined in Proposition 5.2 
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Consider first (8.2). Since T2 = ri then r?(Ti) < < |loge|~5. From (6.3) we see 
that in fact D(a(ri)) ^ 1 <C 100. We now claim that there exists a /xq such that for 
ah Ti < t < Ti + //o we have D{a{i)) < 100. In particular, by (|1.11[) and (2.5) 



D{a{t)) = D{a{n)) + 

< D{a{n)) + r 

J T1 



I • |2 

a 



I • |2 

a 



n 



Ti Jn 



\dtu\^ 
TT |loge| 



< D{a{Ti)) + Cfio^ < 100 
Pi 

for fiQ small enough. 

Next consider (|8.3[). Again r/(ri) < |loge|~5 and so by Lemma 



3.2 



we have 



\\Jiu){n)-Y,dj6a^ 



By Step 2 there exists a /ii > such that for all ti < t < ti + fii 

II J(n)(t) - J(n)(ri)||^_,,, < C(i?,(no))a;(M) < ^ llogep ^so 
for /xi small enough. Furthermore, there exists /X2 such that for ti < t < ri + /i2 

3 



for /Lt2 small enough, where we used (2.5) in the third inequality. Therefore, for 
IJ, = min{;^0) IJ-i, ^2} > and all ri < t < ti + /u we have 



<\\J{u)it)-Jiu)in) 



\w- 



+ \\Jiu){n)-^d, 

^aj{Ti)\\w-l,l 

and D{a{i)) < 100. Hence, we can extend ri to ri + //. 

4. Assume that T > de > ti > T2, then we claim that ti > T2 > Ss- Since r2 > 
we can use (4.8). 

_7_ 1 

|?7| < |loge| 80 I log |loge| 1 2 

then 

_x 1 1 

r/(i) < T/(0) + tC|loge| so |log|loge| |2 < -2?^. 

Hence for all < t < |loge|~6 then rj^t) < ^T)^. In particular, this shows that T2 = ti 

_ 1 

up to 5e < I log e I 6 . From Step 3 we can extend ti and T2 to at least 6e- 
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5. Assume that T > ri > T2 > Ss, then we claim that ti > T2 > T. Note that 
V3 ^ V2 using ( |4.9[ ). Note that 

_x 1 1 

{r]{de))5, <viO) + C6e\loge\ 80 |log|loge| \^ < -P^. 



From Proposition 7.1 we have the differential inequality for the averaged {ri)g^: 

d n^T 11 

-fAri)5^<\{Tj)sAt)\<— sup (r/)^ (s) + |loge| « . 



Using the Gronwall argument from Lemma 8.1 below, we find 



2T 



n 



exp 



for e < eo for aX\ 5e <t <T = min{^^log |loge| , tq}. By (|49]) we find r]{t) < \Ve 
for t in the same time interval. 

□ 

We conclude with the following Gronwall estimate used at the end of the proof of 
Theorem II. 1[ 



Lemma 8.1. Suppose A,B are positive constants and suppose 

d 



dt 



x{t) < A sup x{s) + B 



then 
(8.4) 



xit) < 



se[o,i] 
B 



xq + 



A 



-,At 



Proof. Let m{t) = sup^g^,*] ^{^) then rn{t) = max{i;(t),0} since the maximum can 
increase only if x increases. On the one hand, if rn{t) = x{t) then rh(t) = x{t) < 
Am{t) + B. On the other hand, if m(t) = then rh(t) < Am{t) + B. The estimate 
follows. 

□ 



9. Estimates on the Neumann function and the renormalized energy 

In this section we provide estimates on the Neumann functions that comprise the 
renormalized energy in the Dirichlet case. These estimates will be used both to 



generate long-lived solutions of (1.1) with asymptotically many vortices and to provide 
kernal estimates for the hydrodynamic limit in the next section. The challenge that 
we face is to control how close the aj (t) can get to each other or to the boundary using 
solely information on the renormalized energy W{a). This provides only exponentially 
weak control and leads to dilute concentrations for asymptotically long times. 

The following proposition gives a class of data for which we can be assured that 



the point masses stay separated enough to use Theorem 1.1 on a long enough time 
scale. In particular we need T > ^. 
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Proposi tion 9.1. Fix e and suppose ^^(O) is initial data satisfying the hypotheses of 
Theorem 



1.1 



in the Dirichlet case with ne(0) = +*'^* for ip^, G C and Jq^^* 

III — - III- 
0. Furthermore, assume Pa{o) ^ |log | log |loge| || ^ and n < |log | log |loge| ||4 then 

solution Ue{t) has vortices close to the aj{t)'s (in the sense of Theorem \1-1^ , which 

satisfy 

Pa(t) > I log I log e I 

for all < t < \ log |log e| | ^ . 



' 10 



In order to prove Proposition 9.1 we follow the approach of Sandier-Soret [45j to 
define the renormalized energy in terms of Neumann functions. In particular let 
Nn{x,y) denote the Neumann function which satisfies the following equation 

ANn{-,y) = 6ymn 

(9.1) 1 

duNn{-, y) = drO + -9r9J* on d^l. 
n 

and the limiting Neumann function N{x,y) = Noa{x,y) which satisfies the following 
equation 

AN{-,y) = 5ymn 
^ ■ ^ duN{-,y) = dr9 on d^. 

We also define Hn{x, y) = Nn{x, y) — log |x — y| and H{x, y) = N{x, y) —log |x — ?/| to 
be the harmonic pieces of the Neumann functions Nn{x,y) and N{x,y), respectively. 
Then 

n 

(9.3) W{ai, . . . ,a„) = -vr ^ iV„(aj, a^) - tt'^ Hn{aj,aj), 

j¥=k i=i 

see [3] for example. 

We state the following useful set of estimates: 

Lemma 9.2 (Sandier-Soret [IS]). The Neumann function Nn{x,y) verifies 

(1) Nn{x,y) = Nn{y,x) 

(2) Nn{x, y) = log l^— y| +Hn{x, y) where Hn{x, y) is continuous onQx ilUfi x $7. 

(3) Nn{x, y) = 2 log |x — y| +Hn{x, y) where is continuous on dQ x fiuri x dO,. 
for all I < n < oo. 



In the proof of Lemma 9.2 the authors generate Hn{x, y) in steps. When = Bi 
and duN(-,y) = 1 then H{x,y) = H{x,y), where the complexified H is explicit: 

(9.4) H{x,y) = log\l — xy\. 

For nontrivial d,yNn{-,y) = fn = drO + ^(^9^ one finds Hn satisfies Hf^{x,y) = 
H{x,y) + P{x) + Q{y) where P{x) and Q{y) are harmonic in Bi and bounded and 
continuous up to the boundary. Finally, for simply-connected domain Q let w{z) 
denote the conformal mapping of Q into Bi. Then one finds 

(9.5) Hn{x,y) = Hj^{wix),w{y)) 
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where is defined as fn{w{z)) = fn[z)/w{z). We note again that in our case 

fn{x) = drO + \dr^i.. 

From Lemma 9.2 and the construction of Hn, one can readily prove the following 
estimate on Hn- 

Lemma 9.3. The functions Nn{x,y) and Hn{x,y) satisfy the following estimates 

< C 



(9.6) 



Unix, y) — log [max{|a; — y|, dist(a;, dVt), dist(y, 5r2)}] ^ 

Nn{x,y)<C 

Hn{x,x)<C 

uniformly in n, where C depends on Vt and 93^. 



Using (9.6) we prove the following lemma that provides a lower bound on the 
intervortex distance as a function of the renormalized energy. 

Lemma 9.4. Let W{a) he the renormalized energy for the Dirichlet case. Then 

(9.7) Pa > e-^(-)-^"' 

for some constant C that depends only on Vi and 



Proof. Since the domain is bounded we have from (9.6) 
W{a) = - ^ Nn{aj,ak) - ^ Hn{aj,aj) 

j¥=k j 

> log[min{|oj — Ofcl}]^"'^ + log[min{dist(aj, 50)}]"^ — Cn^ 

> log p-'-Cn^ 

where C depends on if^, and O. Therefore, 



Pa 



>gPa > ^-W(a)-Cn^ 



□ 



This lower bound is sufficient to continue solutions to asymptotically long times. 
In particular we finish with the 



Proof of Proposition 9.1. Let aj{i) be the solution to 



1 



vr 



in the Dirichlet case. Our initial data is chosen to ensure that p^ 
From our assumptions on Pa(Q) and n we have 



log I log e\ 



n 



VF(a(0)) < C < C| log I log |loge| 

Pa(0) 
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Since W{a{t)) < W{a{0)) then from Lemma 9.4 we see that 

5 

< I log |loge| 1 6 

for ah time; hence we can set p^, = \ log | |loge| Finally, this implies a time of 

existence up to 

T = C\ log llog e| I i ^ > C > I log llog .1 I f , 

^ I log I log |loge| II 4 

which certainly satisfies T > K □ 

10. Hydrodynamic limit 

We start by examining the vortex density function uj(t) = lim^^oo ^ 12^=1 ^aj(t) = 
lim„_>oo which is the A^-limit of the cloud of delta functions given by the ODE 
We follow the approach of Schochet flBj and Liu-Xin [HH] for the planar 
Euler equations and Lin-Zhang [34| for the planar time dependent Ginzburg-Landau 
equation. Here we encounter some challenges due to the boundary condition, and this 
leads us to consider interior weak solutions. 

In the following subsection we will relate this ODE limit to the energy density 
^ corresponding PDE limit problem for an appropriately chosen sequence 
e„ with appropriately constructed initial data. 

Proposition 10.1. Consider a sequence of initial data {aj(O)}"^]^, and assume 



for every n. Let aj{t) solve ( 1.11 ) with W{a) in the Dirichlet case. IfuJn = ^ Sj=i 



then resettling time t = nt, we find ujnit) ^{t) in M for all t and uj is a generalized 
interior weak solution to 

d^uj + div (V (Aj^^w) w) = 



with ujQ = lima;„(0) defined by ( 1.18| ). Finally, we have v G Lf^^{Q) where v = 
We first show that the vortex density function conit) satisfies an equation very close 



to ( |L18l ). Following [3], let 

n 

duG = ndrO + 9r(/7jr 
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then G{x) = Yl^=i ^n{x,aj) where Nn{x,aj) is defined by (9.2). If 

n 

Siix) = '^Nn{x,ak) -log|x-aj| 

then from [3] 



fc=i 



-Va,W{a) = VSi{aj) 



where 



V5^(x) = ^VA^„(x,afc) 



n 



VNn{x,aj 



X — ai 



X — a 



"^Nnix, afc) + VHnix, aj). 
Next note that 6a^. = AA^„(-,afc), so for a test function x G C^{Q) and with ujn 



^ E?=i '5a, (t) we have 



\ I 1 " ^ 1 " ^ 



11 11 

-^rr? [Va,^(a)]^ = ^ ^ 5£X(«i)5^5^(ai 

j=l j=i 

1 1 



X] / di,x{x)diSi{x)5a^{x). 



Using the above identity for V5n and 5a {x) yields 



1 /" 11 /■ 

-5t / X^^n(i) = 2 / ^(-Xix) 

J j=i J 



deNn{x, ak) + diHn{x, aj) 



dmdmNnix,aj)dx 



1 

1 1 



+ 

TT n 

^ An + Bn. 



dixix)diNn{x, y)dmdmNn{x, z)ujniy)ujn{z)dydzdx 
dix{x)deHn{x, y)dmdmNn{x, y)ujn{y)dydx 



Following [36] we define the matrix-valued function IC{n,y, z;r]) 

(10.1) ICjkin,y,z,ri) = / r]{x)dx^Nn{x,y)dx^Nn{x, z)dx, 

Jn 
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and after a short calculation, one can rewrite An and Bn as 
1 



(10.2) 



Bn 



27r 
1 

vr 

nvr 



(/Cii - /C22) (n, y, z, (9^^ - d^^) x) ujn{y)i^n{z)dydzdx 



^12 (f^, y, z, dx^dx^x) ujn{y)'^n{z)dydzdx 
dmdm. {deX diHn{x, y)) Nn{x, y)ujn{y)dydx. 



We will show that as n — )• c«, Bn converges to zero and vln's converge to the form of 
the generalized weak solution. However, in order to complete the proof, we prove two 
technical lemmas on the KLjt and the vorticity maximal function (defined below). 

The following lemma provides estimates on the Neumann functions Nn{x,y). In 
particular we have: 



Lemma 10.2. The matrix functions )Cjk{n,y, z,rj) defined in (10.1) safisfies the fol- 
lowing estimates for y,z ^ ft and rj E Cq'{Q): 

(10.3) |(/Cii-/C22)(n,y,z,r/)| < C 

(10.4) |/Ci2(n,y,z,??)| < C 

(10.5) |/Cii(n,y,z,??)| + |/C22(n,y,z,r/)| <21og|y-z| + C 
where C depends only on rj, ip^,, and 0,. Finally, we have the bound 



(10.6) 



V^i/n(x,y) 



< 



C 



dist{y,dny 



where C depends on k, ip^, and 



Proof. These estimates are similar to ones found in Delort [9j and Evans-Miiller |16j 
for the associated Greens function on M^; therefore, we sketch the proof of (10.4) 
following the argument of fT6]. The proofs of (10.3) and (10.5) can be established by 
similar adjustments of arguments in |16j . 

To prove (10.4) one needs to examine the behavior of the gradient of Hn{x,p) = 
Nn{x,p) — log \x—p\ defined via (9.4) and (9.5). Since the test function ry has compact 
support away from the boundary, it follows that dx.Hn{x, •) is bounded for all x on 



the support of r/, as in the proof of Lemma 9.2 
We can now write 

'{x - y)i 



l/C 



12 



< 



r]{x) 



+ dx^Hn{x,y) 



+ 



rj{x) 



\x - y\ 
{x - y)i {x - z)2 
\x — \x — 

dx-,Hn{x,z)j^ 



dx 



+ 

y)2 



r]{x) 



y\ 



+ 



^ + dx2Hn{x,z) 
{x-y)i r^ , ■ 

I 19 '-'X2^n\X , Z) 

\x — 

r]{x) [dx^Hnix, z)dx2Hn{x, z) 



= I^+l2 + h + h 

Using the support of r] and the explicit estimates in the proof of Theorem 1.1 of [16], 
it follows that h < C. I4 < C due to the uniform bounds on VxHn{x, •) for x having 
compact support away from the boundary, where C depends on the distance of the 
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support to the boundary. Finally, we consider the the bound on I2 and /s, which can 
be handled by identical bounds. Due to the uniform bound on VxHn^x, •) away from 
the boundary, we have 



< 



diam(n) 



{x - y)i ^ 

I in C'X2^n\X 1 Z ) 

\x — 
dr < 1. 



1 



supp(?;) I''^ y\ 



Combining the estimates yields (10.4). 



□ 



Define for any Radon measure ^ the maximal vorticity function Mr{n) of DiPerna- 
Majda [lO] 



Mr{n)= sup / \fi{y,t)\dy 
xenfi<t<T J Br-(x)nn 



{x)nn 

for < r < As in [371 EH EH we prove a decay estimate on Mr(w„) below in order 
to pass to the limit in the main term An- In particular we have: 



Lemma 10.3. Suppose {aj{t)Yj^-^ arise from the hypotheses of Proposit 
we can bound 



ton 



10.1 



then 



n 



for all n and all r < 1. Furthermore, 

Mriu;) < 



Proof. Following the structure of the argument in [36] we have for some positive 
integer < n, 

1 2 



logr|M^ iujn{t)) = I logr| 



1 



n 



#{aj{t) G Br{x)nn} 



I kx{kx 1) , 1 1 \kx 
iogr| ^ h I logr| — 



< 



I2 Yl [-Nn{aj,ak) + C] 



<1 + 



\aj-ak\<r 

log r| 



I loE r\ 



n 



n 



MriuJn) 



where we used Lemma 
bound follows. 



9.2 



Lemma 



9.3 



and X^jyfc \ < r? . Since Mr{ujn) ^ 27r the 



For the bound on Mr{uj) we have for x £ where x = 1 on Bj.{x) and x = on 
i^\B2r{x)., X is chosen where (a;)nc'^(*) maximal, then 

Mr{u:) < / = lim / x^n < lim M2r{uJn) < -7=- 

J n-^ooj n^oo V | log 



□ 
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Proof of Proposition 10.1. We now examine the convergence behavior of An and 



From Lemma 10.2 one can fohow the arguments of |461l36j to estabhsh the convergence 



of All. Looking at A\ and taking x £ C'o°(^) and setting rj = (5^^ — S^^) we have 



" 27r 
1 

~ 2^ 



{\y-z\>r}r]n 



(/Cii - /C22) (n, y, z, T]) ujn{y)uJn{z)dydzdx 



{a<\y-z\<r}r\U 



(/Cii - /C22) (n, y, z, 7]) LOn{y)i^n{z)dydzdx. 



Since (/Cn — /C22) y, -z, f]) is continuous in each variable and bounded in the region 
then that term converges to 



1 

'2^ 



{\y~z\>r}nn 



(/Cii — /C22) (00, y, z, rj) uj{y)ijj{z)dydzdx. 



On the other hand in the second region we have 



1 

2^ 

< C 



(/Cii - /C22) y, z, rj) uj{y)uj{z)dydzdx 



{Q<\y-z\<r}r\n 

iOn{y)uJn{z)dydz 



{a<\y-z\<r}r\U 



^ ll'^nllLi(n) 



{|z|<r}nn 



ljJn{z)dz 



and by Lemma 10.3 the term goes to zero as n — )• cx) and r — )• 0. This implies 
A^^ — >■ . The convergence of is much easier since the kernal is continuous on the 
entire domain. 

Next, we show that Bn — 0, and here we crucially use the compact support of the 
our test function x- Bn consists of three terms, depending on where the derivatives 



hit. We consider the worst case in which all derivatives hit Using (10.6) we get 

dixdmdmdiHnix, y)Nnix, y)ujn{y)dydx 



1 

nvr 

< - ||w„||ii 
n 

as n — )• 00. The rest of the terms of Bn are estimated in a similar fashion. 

Finally, we can prove the estimate on the kinetic energy in the fashion of Liu-Xin 
[36] . As in [36] one can use the decay of Mr{uj) — )• to prove that 



(10.7) 



{\y-z\<r}nn 



log \y - z\uj{y)uj{z)dydz < 1. 
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Then for K a compact set in take a nonnegative test function x ^ with 
X = 1 on K. Then 



V < XV 



K 



{|j/-2|<r}nC 



(/Cii + /C22) uj{y)uj{z)dydz 
(/Cii + /C22) uj{y)uj{z)dydz 



+ 

'{\y-z\>T}nn 
A + B. 



(/Cii + /C22) uj{y)io{z)dydz 



Since B is away from the singularity, then we see immediately that B is bounded. 

□ 



The bound on A follows from (10.5) and (10.7). 



We are now in position to establish the hydrodynamic limit. The primary task is 
to approximate the initial data in a suitable way by quantized vortices that satisfy a 



good energy bound. Then we can use Proposition 10.1 



Proof of Theorem 1.3. We first approximate initial data for < wq G A1 H iSf ^(J^) 



in a suitable way so that we can use both Theorem 1.1 and Proposition 10.1 



1. Assume suppcjo C with dist(J7,9f]) > C > 0. We then cover our set O with 
nonoverlapping squares {Qj}, where 

Qj = j'th square of side-length h, 

so there exist 0{h^'^) squares Qj that cover il. We then set 

(10.8) /i = n-3 

so <C n. We now define 

where the Wg ^ are set below. Next, set 



n 



and Uj = [n^J , then | nj — | < 1 and 



Wo 



(10.9) 



< h 



-2 



n2 . 



Since ujq has compact support, then for all h < ho = /lo(^) small enough, if QkCidi^ 7^ 
then rij = 0. If we set n = rij then 

n — Cn^ <n<n 

so n — )• 00 in the same rate as n — )• 00. We can then use n instead of n in the 
discussion below; however, we relabel n as n for simplicity. 
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Next, we slice Qj into rij thin rectangles of equal width. They will be aligned 
vertically and horizontally in alternating sequence, see Figure [2] In the center of 



Qk 



Q 



k+1 



Figure 2. Construction of the rectangles in the QkS 

each of these subrectangles we label points {flg • • • '^ojli ^o the distance between 
neighboring points is — . Finally, we let 



(10.10) 



k=l 



where the are defined above. In the worst-case scenario all vortices are located 

in a single cell with intervortex distance O(^) n~4, and we will need to check that 
this conforms to the correct bound on Pa(o)- 

We claim that Wg — ?■ loq in A^(il). Let fu denote the average of / on U . Then for 

X e C^{^), 1/^ X H - ^o)| < Eq,. Sq^ (x - XQ,)K - ^o) +Eq, \xqj I - wo 

as n — )• oo from ( 10.8 ), ( |10.9 ), and the continuity of x- Therefore, oj^ — )■ ojq 'm.M.[Q). 
2. Finally, we claim that 

(10.11) ^"K^.<fc)<l• 

Since the support of lies in a compact set away from the boundary, then 

mm{dist{alj,dQ)} > C > 
uniformly in n. Hence, we have |-ff„(og •, Og j)| < C uniformly in n. In particular, to 



establish (10.11) it is sufficient to prove 

1 



(10.12) 



< 1. 
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We subdivide the sum into those vortex interactions arising from the same Q/s and 
those that arise from differing QkS, 



= A + B 



We now consider the sum A. Concentrating on a single Qj, assume without loss of 
generality that the subrectangles are vertical and j is located at the origin. Then the 
vortices in this square are located along the x-axis with x values at {0, A, 2A, . . . , {nj — 
1)A}, where A = Summing over the log interactions yields 



- log \aoj - a{j\ = - [{uj - 1) log |A| + (nj - 2) log |2A| + • • • + log |(nj - 1)A| 



< ^logA-i < ^log/i-^ + ^logn 

< 3n2log/i-\ 



since log/i ^ = jlogn. Now summing over the j's yields and using that ^ < 
irlg^^o, we get 



^ = E E i°g \<i - «o,,i ^ ^ E 

j kf^e j 
<Vlog/i"W ujoiy) uJo{z) 



^"E/ / '^o(y)log|y - z|a;o(2:)dydz 



Next we bound B. Let pj denote the center of the square Qj. Due to the alternating 
alignment of the subrectangles in Figure j2| we see that log [oq ^ — Oq fcl ~ \Pj ~ Pk 



< 
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C, even for neighboring squares. Therefore, 



jyfc i=l £=1 



- ~ ^ E E E ~ ^fcl = ~ E ^i"^ li?* - Pfc 



j^k i=l i=l 



< I - ^log\pi - pk\ wo(y) / ujo{z) 



< 



^o{y) log |y - z\u}o{z)dydz 



Combining A and B together we find 



^ + B <l - I I LOQ{y)\og\y - z\u)Q[z)dydz 
Jn Jn 



< 1 + ||t^o|li/-i(n) 



/^supp(i^o)(y) / ^og\y - z\uJo{z) 



< 1 + llwoll^-i < 1, 



where fiq is the characteristic function on Q. 

3. Now we complete the proof of the hydrodynamic limit. Set e„ such that n 
I log I log I log 1 1 1 4 and 



n 7r|loge„J 



Given the initial measure ujq, we build our initial data tte„(0) with vortices at {«oj} 
as generated above. The intervortex distance is no worse than 



h _3 III — 2 

Pa(o) >C->Cn 2 > C log log log e„ » > log log log 



which satisfies the requirements of Proposition 9.1 



We now take a test function x G then setting 



vr loge„ 



it) 



dx 
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we have 



1 

n Jn 
1 



UJn{t)xdx 



Jn 




+ A(t)-A(s) 
1 " 



n 



t n 

t 



Vx(a,(r))-d,-(r) + A(t)- A(s) 

-- r Vx(a,-(r)) • Va^H^(a) + A(t) - A(s) 
vrn 7^ 



1 

47r 
1 

2^ 



((^1 - ^l)^) (^^1 - V2) dxdydzdr 



did2X V1V2 dxdydzdr 



as n — >• 00, using Proposition 10.1 Here we used the simple estimate 



IIVxIIj 



Vrl log Er 



<1 

9 



0. 



w- 



Since this is true for every positive s < t, we complete the proof of Theorem 1.3 



□ 
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